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ABSTRACT 

The  deformation  of  an  unstable  baroclinic  wave  as  a  result  of  mean 
horizontal  wind  shear  is  studied  with  the  use  of  a  five- layer  primitive 
equation  model.   Two  basic  wind  profiles  were  tested:   1)  uniform  zonal 
flow,  2)  jet  type  zonal  flow.   The  jet  profile  produced  a  cold  front 
and  a  warm  front,  with  the  warm  front  being  the  more  intense.   The 
uniform  zonal  flow  produced  a  cold  front.   The  cold  fronts  were  pro- 
duced by  horizontal  shearing  deformation  fields  and  the  warm  front  was 
produced  by  horizontal  stretching  deformation  fields.   The  formation  of 
the  warm  front  was  also  enhanced  by  strong  horizontal  convergence.   The 
occlusion  process  was  not  observed  in  the  experiments. 
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I.   INTRODUCTION 

The  concept  of  atmospheric  fronts  has  been  accepted  for  many  years 
yet  at  this  time  there  is  no  mathematical  theory  which  can  explain  the 
formation  of  a  front.   Williams  [1967]  has  indicated  that  the  basic 
problem  is  finding  a  process  which  can  produce  a  small-scale  motion 
field  from  one  which  initially  contains  only  large-scale  motions. 

A  number  of  studies  have  been  conducted  evaluating  mechanisms 
which  are  contributing  factors  in  the  formation  of  a  front  from  an 
initial  large-scale  baroclinic  wave.   Stretching  deformation  wind 
fields  have  been  studied  by  Stone  [1966]  and  Williams  and  Plotkin 
[1968],  while  shearing  deformation  wind  fields  have  been  studied  by 
Arakawa  [1962]  and  Williams  [1967].   Stone  [1969]  indicates  that  these 
effects  should  be  of  about  the  same  importance.   Arakawa  [1962], 
Edelmann  [1963],  Williams  [1965,  1967]  and  0kland  [1969]  have  studied 
atmospheric  frontogenesis  with  vertical  deformation  fields  with  the 
general  result  that  a  realistic  front  is  produced  within  a  reasonable 
time  period.   Bushby  and  Timpson  [1967]  and  Bushby  [1969]  have  studied 
the  dynamics  of  fronts  and  frontal  rainfall  with  a  primitive  equation 
model  utilizing  real  data  with  good  results  in  a  twenty-four  hour 
forecast  of  rainfall. 

The  purpose  of  this  study  is  to  investigate  the  deformation  of  a 
unstable  baroclinic  wave  as  a  result  of  horizontal  shear  in  the  mean 
wind  field.   The  primitive  equation  model  is  basically  the  model  used 
by  Harrison  [1969]  with  a  slightly  modified  lower  boundary  condition 
and  different  geographical  and  theoretical  areas  of  application. 
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II.   THE  PRIMITIVE  EQUATIONS 

4 

Consider  the  primitive  equations  in  the  following  form,  with 
heating,  friction  and  curvature  terms  neglected  and  with  the  coriolis 
parameter  a  linear  function  of  y: 

^=   -   L(u)    -£  +   fv,  (1) 

_  =   .   L(v)    .  _  .   fu>  (2) 

ff  =   -    L(9),  (3) 

c\jd  _       fdu        cV\  ,, . 

5  --  As  +  w>  (4) 


^p       p  ^p  n?0. 


u*)M>  <*> 


where 


L(s)   -^   (us)  +^   (vs)  +^   (cds), 


and  s   =  u,   v,    0. 

At    p   =   p 
r        ro 

St  "  -  "S  (U0)  -  ~y    (V0)  +  0  VS  +  SyV  "  <°  ^  •      (6) 
Equations  (1)  through  (3)  and  equation  (6)  are  prognostic 
equations.   Equations  (1)  and  (2)  are  the  momentum  equations, 
equation  (3)  is  the  first  law  of  thermodynamics  and  equation  (6) 
is  the  lower  boundary  condition.   The  diagnostic  equations  (4)  and 
(5)  are  the  mass  continuity  and  hydrostatic  equations,  respectively. 

Equation  (6),  the  lower  boundary  condition  was  obtained  by 
setting  d0/dt  =  gw  =  0,  which  implies  that  w  =  0  at  the  lowest  level, 
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Phillips  [1963]  has  shown  that  the  error  resulting  from  the  neglect  of 
vertical  motion  at  the  lowest  level  is  very  small. 
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III.   FINITE  DIFFERENCES 

Equations  (1)  through  (6)  were  solved  numerically  by  introducing 
finite  differences  in  x,  y,  p  and  t.   The  space  finite  differences 
employed  were  a  variation  of  a  scheme  constructed  by  Arakawa  [1966] 
[also  see  Lilly,  1965].   The  basic  scheme  is  designed  to  conserve 
total  energy  and  also  to  conserve  the  squares  of  u,  v  and  9  with 
respect  to  the  advection  terms. 

The  vertical  domain  is  divided  into  five  layers  with  the  arrange- 
ment of  variables  as  shown  in  Fig.  1. 

mb  Variable  Level  (k) 

0    CO  ■  0   10   


200    8  

400    6  

600    4  

800    2  

900 u,v,0,0 1 

1000    u,v, 0,0,05  0  

Fig.  1.  Distribution  of  variables  in  the  vertical. 

The  finite  difference  form  of  the  L(s)  operator  can  be  written 
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L(s)i,j,k =  L(ui+l,j,k  +  "i.j.k'  (si+1,j,k  +  st,j,^ 

i,j,k       Ui-l,j5k       Si3j,k       Si-l,jJkJ 

+  [(vi,J+i,k  +  vi,j,k)  <sij+i.k  +  Sj.k* 

-  <"i,i,k  +  vi,j-i,k»  <sij,k  +  si,j-.,k)]Wl' 


The  linear  space  derivatives  are  of  the  general  form 

|j  =  (s.^.  .  ,  -  s.  .  .  .  )/2Ax. 
dx     i+l,j}k    i-l,j,k' 

The  hydrostatic  equation  (6)  can  be  written  in  finite  difference 

form  following  Lorenz  [1960]  as 

(9i    i    k+  9i    1   k+2} 
i,j,k+2  i,j,k  2 


<WK 


ro  ro 


The  form  of  vertical  finite  differences  as  specified  in  the  L(s) 
operator  and  the  hydrostatic  equation  were  modified  at  level  zero  in 
the  following  manner. 

1.  9,  u  and  v  were  linearly  extrapolated  from  levels  one  and 
three  to  level  zero. 

2.  d0/dp  at  level  zero  was  approximated  by  a  one-sided 
difference . 

Centered  time  differences  were  used  throughout  this  study  with  the 
exception  of  the  first  time  step  which  is  forward.   To  prevent  sepa- 
ration of  solutions,  the  program  was  restarted  every  twelve  hours  with 
a  forward  time  step. 


15 


The  following  basic  steps  were  used  in  integrating  the  given 
system  of  equations. 

1.  The  potential  temperature  and  wind  fields  at  the  odd- 
numbered  levels  and  the  geopotential  of  the  1000  mb  surface  are 
prescribed  initially. 

2.  The  mass  continuity  equation  is  integrated  in  space  from 
the  highest  level  to  the  surface  to  obtain  cjd. 

3.  The  hydrostatic  equation  is  integrated  in  space  from  the 
surface  to  the  highest  level  to  obtain  0. 

4.  The  momentum  equations  are  integrated  in  time  to  obtain 
u  and  v  at  time  t  +  At . 

5.  The  first  law  of  thermodynamics  is  integrated  in  time  to 
obtain  0  at  time  t  +  At . 

6.  The  equation  for  the  lower  boundary  condition  is  integrated 
in  time  to  obtain  the  1000  mb  geopotential  at  time  t  +  At . 

7.  Steps  2  through  6  are  repeated  for  the  number  of  time  steps 
desired. 
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IV.      BOUNDARY  CONDITIONS 

The  boundary  conditions  are  specified  in  a  manner  that  conserves 
mass  and  energy  within  the  grid  domain.   At  the  northern  and  southern 
boundaries  the  normal  component  of  velocity  is  zero  and  periodicity  is 
enforced  at  the  eastern  and  western  boundaries.   Fig.  2  illustrates  the 
horizontal  grid  domain. 

i=l  i=M 

j=N  N.  boundary  j=N 


*  i  #  •     •     • 
i 


-  •  •!• 


W.  boundary  E.  boundary 


»  *    * 


*         *>         »        «         *         w  *         «.         V         *        •         * 

i=l  S.  boundary  i=M 

J-l  j=l 

Fig.  2.   Horizontal  grid 
The  east-west  boundary  condition  forces  cyclic  continuity  on  all 
parameters.   This  is  accomplished  by  requiring  the  values  in  the  two 
columns  on  the  right  to  be  the  same  as  the  values  in  the  two  columns 
on  the  left.   This  condition  is  specified  by 

Sl,j  =  VlJ  S2,j  =  SM,j' 

where 

s  =  u,v,0,9. 

The  north-south  boundary  condition  allows  flow  along  the  boundary 

but  does  not  allow  flux  across  the  boundary.   This  is  accomplished  by 

imposing  the  following  condition  on  the  v  term  of  the  advection  operator, 

v#1+v=0  v.    ..  +  v.    „  -   -  0 

i,l  i,2  i,N  i,N-l 


or 


Vi,l   =   "   Vi,2  Vi5N=   "   Vi,N-l 
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The  values  of  geopotential  and  potential  temperature  in  the  outer  row 
are  set  equal  to  the  value  in  the  next  interior  row.   This  is  speci- 
fied by 

s  •  i  =  s  •  i  s.  „  =  s.  „  .       (7 ) 

i,l    i,2  i,N    i,N-l 

i  ■ 

where 

s  -  0,9- 
As  a  result  of  v  being  set  equal  to  zero  between  the  two  outer  grid 
points  along  the  northern  and  southern  boundaries,  values  of  u  and  co 
are  not  required  outside  the  boundaries. 

The  upper  boundary  condition  of  U)  =  0  eliminates  flux  through  the 
top  of  the  domain.   The  lower  boundary  condition  has  been  specified  in 
equation  (6) . 

Various  boundary  conditions  were  tested  (see  Appendix  A)  before 
acceptance  of  the  above  method. 
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V.   INITIAL  CONDITIONS 

The  initial  fields  of  u,  v  and  9  were  specified  in  order  to 
generate  a  baroclinic  zonal  current.   A  small  amplitude  disturbance 
in  the  geopotential  field;  which  is  independent  of  pressure,  was 
superimposed  on  the  basic  zonal  flow.   The  mean  zonal  current  (u)  has 
the  following  forms  ; 

1)  jet    profile 

u  =  5  (p  "  p0}  sin    kwV 

2)  uniform  zonal  profile 

U  "  5J  (P  "  Po) 

A 

where  cki/dp  is  constant.   The  initial  data  fields  were  generated  in 

the  following  steps. 

1.  The  geopotential  at  p  =  p  was  calculated  from 

2.  The  potential  temperature  fields  were  specified  from 

P  xH 


9         r      ap  J  +   o     aP  ' 

A 

where    d9/dp   is   a   constant. 

For   uniform  zonal   flow 
J=foy+goy2/2. 
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For   jet   flow 


J  =   f  y  +  0y2/2 
o 


'<  Cf  -  h  -  er): 


o 
2 


r-     2 


v_  _  w^  cos  rimfN .  vw  sin  /-2nx 


L2         .2  V  w  y       2n  vwyJ 

4tt 


3.  The  hydrostatic  equation  was  then  solved  for  the  geopotential 
at  all  required  levels. 

4.  The  initial  wind  fields  were  then  calculated  from  the  geo- 
potential fields  using  the  geostrophic  approximation.   Thus,  the 
initial  fields  were  in  geostrophic  balance  at  time  zero. 

The  basic  current  is  zero  at  the  lowest  level  and  increases  at  a 
rate  of  5  mps/100  mb  to  a  maximum  of  45  mps  at  100  mb .   Therefore 
du/c3z  decreases  with  increasing  altitude  which  reduces  the  baroclinic 
instability  in  the  upper  layers.   This  approximates  the  effects  of  a 
tropopause  on  the  dynamics . 

The  initial  disturbance  amplitude  was  10  meters  or  .125  mb .   The 
small  initial  amplitude  was  chosen  to  prevent  the  excitation  of  large 
amplitude  gravity  waves  that  occur  since  the  initial  wind  field  was 
geostrophic . 

The  value  of  d9/cip  used  was  -.376  C/mb.   This  value  is  approxi- 
mately one-half  the  value  determined  by  Gates  [1960]  for  winter  over 
the  United  States,  but  it  has  the  advantage  of  forcing  the  fronts  to 
develop  more  rapidly  and  thus  to  reduce  the  computer  time  required. 
The  value  used  approximated  the  standard  atmosphere  quite  well  in  the 
lower  levels,  see  Table  I,  but  is  less  stable  than  the  standard  atmos- 
phere in  the  upper  levels.   The  stratosphere  is  clearly  not  represented 
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by  this  temperature  field,  however  the  dynamical  presence  of  the 
stratosphere  is  represented  by  u. 


Pressure  (mb) 

Height  (meters) 

Jet 

Standard 
Atmosphere 

100 

15325 

16210 

300 

9112 

9160 

500 

5576 

5570 

700 

3015 

3010 

900 

988 

990 

1000 

110 

110 

Table  I.   Average  initial  height  of  constant  pressure  surface. 

The  beta  plane  approximation  was  made,  where  the  coriolis  parameter 
is  a  linear  function  of  y.   The  quantities  f  and  0  were  computed  at 
45  north  latitude. 

Various  combinations  of  grid  size,  time  step  and  wind  profile  were 
tested.   Table  II  is  a  summation  of  these  various  experiments. 
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VI.      RESULTS 

A  A 

All   experiments   were   conducted   using   the   values    of   cHi/dp,    ^0/dp 

and  initial  disturbance  amplitude  as  specified  by  equation  (8)  in  the 
initial  condition  section  of  this  paper.  Table  II  gives  a  listing  of 
the  specific  variations  in  each  of  the  following  experiments. 

Experiment  No.  1.   This  is  the  basic  jet  flow  experiment  with  a 
grid  size  of  200  km  in  x  and  y.   The  initial  1000  mb  height  analysis 
is  given  in  Fig.  3  and  the  initial  height  and  potential  temperature 
analysis  for  the  900  mb  and  700  mb  surfaces  are  given  in  Figs.  4  and 
5.   The  growth  of  the  systems  from  day  one  through  day  four  at  the  900 
mb  surface  is  given  in  Figs.  6  through  9.   By  the  end  of  three  days  a 
cold  front  and  a  warm  front  have  formed  and  the  temperature  gradients 
are  large  by  the  fourth  day.   Note  that  the  warm  front  contains  larger 
gradients  than  the  cold  front.   Large  convergence  (not  shown)  is  associ- 
ated with  the  warm  front.   The  700  mb  analysis  of  height  and  potential 
temperature  for  days  three  and  four  is  given  in  Figs.  10  and  11.   The 
cold  and  warm  fronts  have  the  proper  tilt  in  the  vertical,  with  a 
slope  between  1/125  and  1/250,  which  is  computed  at  day  three  and  for 
both  fronts.   The  growth  of  the  surface  pressure  centers,  divergence 
and  vorticity  at  900  mb  is  given  in  Table  III.   It  should  be  noted  in 
Table  III  the  asymmetry  between  the  maximum  and  minimum  values  of  each 
quantity. 

Experiment  No.  2.   This  is  the  basic  uniform  zonal  flow  experiment 
with  a  grid  size  of  200  km  in  x  and  y.   The  initial  1000  mb  height 
analysis  is  given  by  Fig.  3  and  the  initial  height  and  potential 
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DAY 

1000  mb 
centers  (m) 

Divergence 
(x    10°6) 

Vorticity 

rx  in-5} 

Low 

High 

Max     Min 

Max 

Min 

Initial 

100 

120 

.3      -.3 

.7 

-.7 

1 

90    130 

2,2     -2.4 

1.0 

-1.0 

2 

40    160 

3.7     -7.3 

2.9 

-1.8 

3 

-60    190 

6.5    -17.1 

9.4 

-3.0 

4 

-200    270 

20.7   1  -41.9 

18.9 

-6.6 

Table  III.   Experiment  No.  1. 
temperature  fields  for  the  900  mb  and  700  mb  surfaces  are  given  in 
Figs.  12  and  13.   The  growth  of  the  systems  from  day  one  through  day 
four  at  the  900  mb  surface  is  given  in  Figs .  14  through  17 .   By  day 
three  a  cold  front  has  formed  to  the  west  of  an  apparent  warm  front. 
However  there  is  a  single  vorticity  maximum  in  the  region  and  there 
is  no  evidence  of  a  warm  front  at  higher  levels.   The  front  appears 
,   to  be  an  occlusion  but  it  did  not  evolve  from  two  separate  fronts  with 
associated  vorticity  maxima.   It  was  also  noted  that  in  the  northern 
region  of  the  low  there  was  a  packing  of  potential  temperature  lines 
which  somewhat  resembles  the  warm  front  in  the  jet  flow  experiment. 
The  700  mb  pressure  analysis  of  height  and  potential  temperature  for 
days  three  and  four  are  given  in  Figs.  18  and  19.   The  front  tilts  to 
the  west  in  the  vertical  with  a  slope  between  1/60  and  1/120  at  day 
three.   The  growth  of  the  surface  pressure  centers,  divergence  and 
vorticity  at  900  mb  is  given  in  Table  IV„ 

Experiment  No.  3.  This  experiment  is  the  same  as  experiment  1 
except  that  the  grid  length  was  reduced  to  100  km  in  x  and  y.  The 
results  are  basically  the  same  as  in  experiment  1,  with  the  exception 
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DAY 

1000  mb 
centers  (m) 

Divergence 
(x  10~6) 

Vorticity 
(x  10"5) 

Low 

High 

Max 

Min 

Max 

Min 

Initial 

100 

120 

.3 

-.3 

.7 

-.7 

1 

80 

140 

2.4 

-3.8 

1.0 

-.8 

2 

30 

170 

8.3 

-8.8 

1.8 

-1.5 

3 

-120 

240 

10.9 

-22.0 

10.0 

-2.9 

4 

-550 

340 

23.3 

-43.7 

25.7 

-4.5 

Table  IV.   Experiment  No.  2. 
that  the  growth  of  the  systems  is  now  slightly  faster  with  a  slightly 
larger  vorticity.   The  900  mb  height  and  potential  temperature  analysis 
for  days  three  and  four  are  given  in  Figs.  20  and  21.   Note  that  the 
temperature  gradients  are  larger  in  this  experiment  at  four  days  than 
in  experiment  1.   This  indicates  that  the  larger  truncation  error  with 
the  larger  grid  size  inhibits  the  energy  cascade  to  small  scales. 

Experiment  No.  4.   This  experiment  is  the  same  as  experiment  2 
except  that  the  grid  size  was  reduced  to  100  km  in  x  and  y.   A  front 
similar  to  the  one  in  experiment  2  developed,  but  in  addition,  a 
gravity  wave  with  a  wavelength  of  approximately  188  km  was  super- 
imposed on  the  height  analysis.   The  gravity  wave  was  initially 
observed  at  day  two  close  to  the  north  and  south  boundaries,  and  by 
day  three  had  completely  covered  the  domain.   It  was  apparently  excited 
by  the  large  shear  at  the  boundary  as  generated  by  equation  (7). 

Experiment  No.  5.  This  experiment  is  the  same  as  experiment  1 
except  that  Ay  was  reduced  to  100  km  and  W  was  reduced  to  2000  km. 
Cold  and  warm  fronts  developed  within  three  days  that  were  similar  to 
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the  ones  generated  in  experiment  1.   Superimposed  on  the  long  wave 
pattern  was  a  very  small  amplitude  gravity  wave  with  a  wavelength  of 
approximately  188  km. 

Experiment  No.  6.   This  experiment  is  the  same  as  experiment  2 
except  that  Ay  was  reduced  to  100  km  and  W  was  reduced  to  2000  km. 
The  results  are  the  same  as  in  experiment  2,  but  in  addition,  a 
gravity  wave  of  small  amplitude  and  approximate  wavelength  of  188  km 
was  superimposed  on  the  long  wave  pattern. 

Experiment  No.  7.   This  experiment  is  the  same  as  experiment  1 
except  that  Ay  was  increased  to  400  km  and  W  was  increased  to  8000  km. 
See  Figs.  22  and  23  for  day  three  900  mb  height  and  potential  tempera- 
ture analysis.   By  day  three  a  cold  and  warm  front,  similar  to  the 
results  of  experiment  1,  had  formed  with  the  disturbance  being  limited 
to  the  central  region  of  the  domain. 

Experiment  No.  8.   This  experiment  is  the  same  as  experiment  2 
except  that  Ay  was  increased  to  400  km  and  W  was  increased  to  8000  km. 
Gravity  wave  propagation,  with  an  approximate  wavelength  of  400  km, 
from  the  north  and  south  boundaries  was  too  strong  to  permit  an 
evaluation  of  the  long  wave  system. 

In  none  of  the  experiments  was  there  an  indication  of  a  front  at 
500  mb  or  above.   All  energy  calculations  (see  Appendix  B)  indicated 
that  the  potential  energy  was  three  orders  of  magnitude  greater  than 
the  kinetic  energy  and  that  there  was  a  small  energy  sink  in  the 
system  up  to  four  days . 
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12  J. 


Minimum  contour  100  meters,  contour  interval 
10  meters. 


Fig.  3.   Initial  1000  mb  height  analysis 
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Minimum  contour  800  meters,  interval  30  meters 


Pig.  22.   Experiment  No.  7,  72  Hour  900  mb  height  analysis 
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Minimum  Isotherm  270  A,  interval  5  C. 

Fig.  23.   Experiment  No.  7,  72  Hour  900  mb  potential  temperature 
analysis . 
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VII.   CONCLUSIONS 

The  time  evolution  of  a  baroclinicly  unstable  disturbance  was 
studied  for  two  basic  zonal  currents.   With  horizontal  wind  shear, 
cold  and  warm  fronts  developed  and  are  similar  to  the  results  of 
Edelmann  [1963]  and  0kland  [1969].   With  no  horizontal  wind  shear, 
a  cold  front  formed  which  was  similar  to  the  results  of  Williams 
[1969].   This  cold  front  resembled  an  occluded  front  at  the  surface 
but  did  not  have  the  life  history  of  an  occlusion.   In  general,  the 
warm  -front  formed  from  the  stretching  deformation  discussed  by  Stone 
[1966]  while  the  cold  fronts  formed  from  the  shearing  deformation 
discussed  by  Williams  [1967].   In  the  jet  case,  the  warm  front  was 
further  enhanced  by  strong  convergence  through  the  action  of  the 
vertical  deformation  field.   The  lack  of  frontogenesis  at  upper  levels 
is  explained  by  the  reduction  of  the  vertical  wind  shear  at  those 
levels.   The  classical  occlusion  process  was  not  observed  in  any  of 
the  experiments. 

Further  studies  with  this  model  should  be  conducted  with  the 
addition  of  vertical  and  horizontal  diffusion  and  with  the  addition 
of  latent  heat  of  condensation.   Other  basic  currents  should  be 
tested,  including  those  with  an  explicit  tropopause . 
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APPENDIX  A 

EXPERIMENTS  WITH  BOUNDARY  CONDITIONS 
Before  acceptance  of  the  north-south  boundary  conditions  as 
specified  in  Chapter  V,  various  other  schemes  were  tested  with 
limited  success.   A  sample  of  these  schemes  with  the  results  follows. 

1.  The  boundary  was  placed  at  the  outer  grid  points  with  the 
values  of  all  variables  determined  from  a  linear  extrapolation  from 
the  next  two  interior  grid  points.   The  integration  became  unstable 
after  approximately  120  six -minute  time  steps. 

2.  The  boundary  was  again  placed  at  the  outer  grid  points  and  a 
special  one-sided  difference  modification  was  made  to  the  basic  diag- 
nostic and  prognostic  equations.   Strong  unrealistic  wind  shear 
developed  along  the  boundary  and  the  numerical  integrations  eventually 
became  unstable. 

3.  Cyclic  boundary  conditions  were  applied  in  the  north-south 
as  well  as  in  the  east-west  directions.   The  discontinuity  in  the 
geopotential  field  forced  the  use  of  a  complete  wave  in  both  the 
north-south  and  east-west  directions.   This  dcvbled  the  number  of  grid 
points  in  the  domain  and  made  the  size  of  the  computer  program  too 
large  for  practical  use. 
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APPENDIX  B 

ENERGY  COMPUTATIONS 

In  pressure  coordinates,  the  total  potential  energy  may  be 

written  as 

Po 

EP  =  -  J     C   T  dp  dA,  B.l 

'  A  o 

where  "dA"  represents  an  area  increment.   Now,  substituting  the 
definition  of  potential  temperature  in  (B.l),  the  expression  for 
the  total  potential  energy  becomes, 


EP  =  lJJ   vQr)  dpdA-  B-2 

A  o         ° 
In  pressure  coordinates,  the  total  kinetic  energy  is  written 


as 

Po 


=  y~    \     f  (u2  +  v2)  dp  dA.  B.3 


EK 

A  o 
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13.  ABSTRACT 


The  deformation  of  an  unstable  baroclinic  wave  as  a  result  of  mean 
horizontal  wind  shear  is  studied  with  the  use  of  a  five-layer  primitive 
equation  model.   Two  basic  wind  profiles  were  tested:   1)  uniform  zonal 
flow,  2)  jet  type  zonal  flow.   The  jet  profile  produced  a  cold  front  and  a 
warm  front,  with  the  warm  front  being  the  more  intense.   The  uniform  zonal 
flow  produced  a  cold  front.   The  cold  fronts  were  produced  by  horizontal 
shearing  deformation  fields  and  the  warm  front  was  produced  by  horizontal 
stretching  deformation  fields.   The  formation  of  the  warm  front  was  also 
enhanced  by  strong  horizontal  convergence.   The  occlusion  process  was  not 
observed  in  the  experiments. 
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